Velocity weakening and possibility of aftershocks in nanofriction experiments 
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We study the frictional behavior of smaU contacts as those reahzed in the atomic force microscope 
and other experimental setups, in the framework of generahzed Prandtl-Tomhnson models. Partic- 
ular attention is paid to mechanisms that generate velocity weakening, namely a decreasing average 
friction force with the relative sliding velocity. The mechanisms studied model the possibility of 
viscous relaxation, or aging effects in the contact. It is found that, in addition to producing veloc- 
ity weakening, these mechanisms can also produce aftershocks at sufficiently low sliding velocities. 
This provides a remarkable analogy at the microscale, of friction properties at the macroscale, where 
aftershocks and velocity weakening are two fundamental features of seismic phenomena. 

PACS numbers: 



I. INTRODUCTION 

In the last years, experimental achievements have 
given the possibility to study friction at the microscale 
(nanofriction) , which has renewed the interest in a variety 
of theoretical and applied aspects of tribology. The pro- 
totypical experimental setup is that of the atomic force 
microscope (AFM) where a nanometer size tip can be 
brought into contact with a substrate, and the frictional 
behavior of the system can be studied in detail^. Other 
related instruments, like the surface force apparatus^ 
have expanded the possibilities to cover other regions 
of parameters, allowing for instance to study somewhat 
larger contact areas, in the range of microns. 

The properties of nanofriction differ substantially from 
those at a macroscopic scale. One reason for this is that 
macroscopic friction originates typically in a set of sparse 
contacts located at different spatial positions, and the ob- 
served behavior is a non-trivial combination of these in- 
dividual contacts. The now available experimental tech- 
niques give the opportunity to study one such individual 
contact in detail. 

One not fully attained goal in the field of friction is pre- 
cisely to correlate the results of well controlled nanofric- 
tion experiments with those of friction at a macroscopic 
scale. For instance, in the largest scale realization avail- 
able, namely that corresponding to seismic phenomena, 
friction displays features that seem to be typical of this 
large scale. One of them is velocity weakening (VW). 
This refers to a friction process in which the average fric- 
tion force is a decreasing function of the relative velocity 
of the sliding bodies. It is widely accepted that VW is 
a necessary condition for seismic phenomena to look the 
way they do^. Although on a phenomenological perspec- 
tive velocity weakening is successfully modeled by the 
so-called rate and state equations'*'^, there is no general 
agreement on a consistent first principle description of 
the phenomenon. It is however agreed upon the fact that 
VW must be induced by some kind of mechanism that 
tend to strengthen the contact between the two sliding 
bodies if they remain in the same relative position for a 
long time. 



Probably the most typical feature of friction at the 
seismic scale is the existence of aftershocks (ASs). They 
occur as secondary shocks after some large earthquake. 
ASs occur exclusively due to an internal dynamics of the 
sliding bodies, and are not directly related to the driving 
mechanism. ASs have been shown to be compatible with 
rate and state equations displaying VW— . In addition, in 
recently studied numerical models of seismicity, internal 
relaxation on the plates have shown to give rise at the 
same time to VW and ASsjii^ 

The dependence of friction force on the relative veloc- 
ity between the two sliding objects is routinely measured 
in nanofriction experiments. In many cases, this depen- 
dence follows an increasing dependence on velocity. In 
a minority of cases, VW is observed. VW is well estab- 
lished at the scale probed by the surface force apparatus, 
where it has been directly measured^J^. Single contact 
AFM measurements have sometimes also found a VW 
friction law. 

The physical mechanisms that have been proposed as 
generating VW in nanofriction experiments typically in- 
clude the existence of some internal dynamics in the con- 
tact that strengthens the tip-substrate interaction as a 
function of time. For instance, inii it was argued that in 
the AFM, capillary condensation of water molecules oc- 
cur at the contact between tip and substrate. The VW 
effect observed in this case is suggested to be due to the 
reduced ability of the system to maintain a stable capil- 
lary neck as the sliding velocity increases. On the other 
hand, ini^, VW was associated with surfaces capable of 
forming hydrogen bond networks, in a process that pro- 
ceeds according to some internal time scale. Recently, 
Carpick et al.^^ have found evidence that micro-contacts 
between solids can become chemically stronger as a func- 
tion of time, an effect that can potentially produce VW. 
In any case, VW is a manifestation of the existence of 
some internal dynamics of the contact, with its own time 
scale. This time scale will ultimately determine the ve- 
locity range for which VW is observed. 

To our knowledge ASs have not yet been observed in 
any nanofriction experiment. Actually, the only context 
where ASs are well characterized is that of seismic phe- 
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nomena. Observation of ASs requires first of all of a non- 
smooth sliding behavior, in such a way that relative dis- 
placement of the sliding materials occur discontinuously, 
and a sufficiently low relative velocity in order to clearly 
separate the time scales of main shock (produced directly 
by the driving) and ASs. In the seismic context the dis- 
continuous events are the earthquakes. In nanofriction 
these are the 'slip' stages of a stick-slip processjiA Yet, 
ASs are a very special type of abrupt sliding events. They 
are events that are caused by previous slip events, and 
that are unrelated directly to the driving. This means for 
instance, that if driving is stopped at some time to, there 
can still be some events occurring at t > to which are af- 
tershocks of those that occurred at t < to. Note that an 
unambiguous classification of an event as an aftershock 
requires a sufficiently slow sliding that does not trigger 
additional events in the mean time between the original 
event and its ASs. Also, it is clear from this description 
that ASs require the existence of some internal dynamics 
of the contact, in addition to the dynamics provided by 
the external driving. 

The simplest model that describes a non-uniform, 
stick-slip dynamics in friction experiments is the Prandtl- 
Tomlinson (PT) mode l^^'^^ , which we briefly review in 
the next section. In its simplest form, the PT model 
does not display VW. In this paper we analyze two vari- 
ants of the basic PT model that produce VW, and also 
ASs. They describe two different implementations of re- 
laxation mechanisms in the system. One can be applied 
(although not exclusively) to describe a viscous damping 
in the tip of the AFM, and the second models the possi- 
bility of aging effects in the contact region of the sliding 
objects. We show clearly how ASs become progressively 
more defined as the relative sliding velocity decreases. In 
addition to the constant velocity setup, we also study the 
response to sudden changes in the driving velocity, look- 
ing for typical signatures in the transient behavior of the 
friction force such as a stress peak at sliding initiation, 
and a gradual relaxation of stress upon sliding stop. 

The paper is organized as follows. In the next sec- 
tion we briefly review the basic phenomenology of the 
Prandtl-Tomlinson model that is relevant to our analy- 
sis, and discuss in detail the motivation of the two vari- 
ants we will study. In Section III we show the results 
of the numerical simulations in a constant velocity set 
up, emphasizing the similarities and differences between 
the two, and showing in particular how ASs appear. In 
Section IV we give the results for some non-constant ve- 
locity protocols. Finally, In section V, we summarize, 
and conclude. 



II. THE PRANDTL-TOMLINSON MODEL, AND 
SOME VARIANTS 

A basic understanding of the origin of friction in 
nanocontacts is provided by the Prandtl-Tomlinson (PT) 
modeli^ii^ in which the contact is represented by a sin- 



gle variable x evolving under the effect of some interac- 
tion potential W{x), that can be written as the prod- 
uct of a nominal distance between the two surfaces h{x) 
and the normal force applied to the contact Fn, namely 
W{x) = h{x)FN. The variable x is externally driven 
through an elastic coupling by a prescribed evolution 
U{t). The PT model adapts particularly well to the 
case of the AFM tip-substrate interaction, where it is 
usually considered that W{x) reflects the atomic corru- 
gation potential. However, It has been used also to de- 
scribe more general situations by using an appropriate 
disordered W{x) function, and having in mind that the 
distance between typical minima of the W{x) and the 
overall values of corrugation h{x) are now defined by the 
roughness associated to the substrate. 
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FIG. 1: The basic configuration of the Prandtl-Tomlinson 
model. The contact, represented by a single horizontal co- 
ordinate X, is pushed down with a force Fn over a rough 
substrate defined by h{x), and driven through a spring by an 
externally applied U{t). 

In Fig. [T] we sketch the basic configuration of the PT 
model. The evolution equation for x{t) consists of a de- 
terministic part associated to the forces applied to the 
tip, and a stochastic part originated in the integrated ef- 
fect of many other "rapid" degrees of freedom that are 
not explicitly taken into account in our description, which 
concentrates only in the "slow" variable x. The effect of 
stochastic forces is twofold. On one hand they produce an 
additional term on the dynamic equation of the form of a 
white external noise with an intensity that is a function 
of temperature. On the other hand they are responsible 
for the appearance of "friction" terms, namely additional 
forces that are proportional to minus the time derivative 
of the slow variables in the problem. By introducing 
these terms, we arrive at the basic equation describing 
the Prandtl-Tomlinson modelii: 



Tjox ={U- x)ko + V2%T0(t) (1) 

here U (t) represents the driving, which for a constant ve- 
locity situation is U{t) ~ Vt, and (j){t) is the stochastic 
force associated to the fluctuating environment. The in- 
stantaneous friction between the surfaces can be consid- 
ered to be measured by the stretching of the ko spring, 
and the average friction force can be evaluated as the 
temporal average of this quantity, namely 

Ff^ = ko{U-x) (2) 



We assume that the substrate potential term has some 
typical distance between minima a, and typical height ho. 
We now give a brief overview of the phenomenology as- 
sociated to this equation, and then discuss some possible 
variants that generate VW and other associated effects. 

We first describe the zero temperature case. The be- 
havior of the system depends crucially on the ratio be- 
tween fco and cPW/dx"^. For fco > d'^W/dx'^\ijiQ_x — 
h-oFj^/a^, there is a single equilibrium position of the 
tip for any value U of the driving. This means that the 
dynamics is quasi-static for a vanishing driving velocity, 
and that the friction force vanishes in this limit. In the 
opposite case fco < '^^^/'^^^ Imax there are metastable 
configurations of the tip, and abrupt transitions among 
them, as driving proceeds. This produces a finite friction 
force in the limit of driving velocity V ^ 0. From now 
on, this is the case we consider. 

At large driving velocities, the friction force becomes 
proportional to velocity, with a pre-factor given by the 
friction coefficient 770- A crossover velocity that is or- 
der rjQ^ can thus be defined, separating the high veloc- 
ity regime of friction proportional to velocity, and a low 
velocity regime of rather constant friction force. We can 
formally describe this low velocity limit by taking 770 — >■ 
in Eq. 

Temperature produces important changes in the veloc- 
ity dependence of friction at very low velocities. In fact, 
at finite temperature, the system has the possibility to 
overcome the energy barriers before the pulling takes the 
system at the mechanical instability point. Since close 
to these instability points the energy barriers are much 
smaller for forward than for backward movements, the 
thermal jumps always reduce the stretching of the fco 
spring, reducing the friction force. The effect becomes 
more important as the velocity is reduced, leading typi- 
cally to a logarithmic increase of friction force with ve- 
locity. A detailed account of this behavior can be found 
for instance ini^. 

Temperature effects are typically important in AFM 
measurements, at least at the lowest velocities, because 
of the small size and mass of the tip. They are less im- 
portant when the size of the contact is larger, as in the 
surface force apparatus. Here, in order to describe in 
the simplest possible form the effect of terms generating 
VW, we will work in the T = limit. Also, as we are 
mostly interested in the limit of very small velocities, we 
will consider from now on the limit 770 of Eq. [T] 

We now introduce two different modifications to the 
Prandtl-Tomlinson model that produce VW. The first 
mechanism replaces the single spring fco driving the sys- 
tem, by a viscoelastic element as displayed in Fig. [2][a). 
This produces a viscous damping when the system is 
driven, and introduces an additional degree of freedom, 
namely the variable 9, and an internal time constant in 
the problem r = (fci -I- fco)77. In the particular case of the 
AFM setup, the present mechanism corresponds to in- 
corporate viscous damping effects in the microscope tip. 
This is one example of the generalizations of the Prandtl- 
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FIG. 2: Two variations on the basic configuration of Fig. 
[1] that produce VW and ASs. In (a), the driving spring is 
complemented by a viscoelastic element. In (b) the tip is 
formally divided in two parts, each of them interacting with an 
independent rough substrate. The two parts are connected by 
a viscoelastic element (represented here by the thicker line). 
See the text for the equations defining both models. 

Tomlinson model known as the "two masses-two springs" 
model^^. The full phenomenology of the two masses- two 
springs model is quite complicated, particularly if mass 
terms for the tip and cantilever are included. Here we 
will not introduce mass terms, and the equations that 
describe this variant (that we call model A) read 

7702; = ~ ^ ~ ~ 
77^ = {U - X - e)ko - 0ki (3) 

A similar kind of dissipation mechanism in the tip for 
AFM experiments has been discussed by Reimann and 
EvstigneexiSSiSi, who have shown how it can produce a 
VW effect. However, they did not present a detailed 
description of the sliding process, and in particular they 
do not refer to the possibility of ASs due to the viscous 
damping. 

Model A describes a modification of the FT model in 
which the new internal dynamics, is introduced "in se- 
ries" with the contact, which itself is not modified. A 
second possibility is to introduce relaxation right at the 
contact. In order to do this, it seems necessary to go 
beyond the description of the contact in terms of a single 
slow variable x. In this respect we note that Evstigneev et 

22,23 ]^QyQ provided evidence that in some cases the re- 
sults of friction AFM experiments are incompatible with 
an interpretation in terms of a single variable characteriz- 
ing the evolution. They have suggested that aging effects 
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have to be included, and introduced a model with an ad- 
ditional degree of freedom describing a "reinforcing" of 
the contact if this has lasted more than a certain time. 

Here we will make an effective description considering 
that the tip-substrate interaction region is not point like, 
but has some spatial extension. This will allow different 
parts of the contact region to accommodate relatively to 
each other to generate a stronger contact, if enough time 
is available. A spatially extended tip-substrate contact 
forces us to consider several variables describing the in- 
teraction in different parts of this region. Trying to keep 
the model in the simplest non-trivial form, we will con- 
sider only two parts of the contact region, described by 
two variables xi and X2. They will respectively move 
on disordered potentials Vi and V2 (producing forces /i 
and /2) which for simplicity will be assumed to be un- 
correlated. The average position of the tip is defined as 
X = {xi + X2)/'2. The relaxation effects will be modeled 
through a viscoelastic response of the variable 9 = x\—xi 
that defines the internal state of the contact. A picto- 
rial sketch of the model, called model B, is shown in Fig. 
[21(b). The evolution equations for x and 6*, that com- 
pletely define the dynamics of this model are as follows 

ryoi = (C/-a;)fco-f /l(x + ^?/2) + /2(a;-^^/2) 
(r/ + r/o)^ = /i(a; + 0/2)-/2(x-(?/2)-fci0 (4) 

When applied to AFM experiments, where there is an 
obvious asymmetry between the two sliding objects, the 
presentation we made of model B suggests that we are 
introducing a relaxation mechanism in the tip, and that 
substrate is inert. This is not necessarily so, as alterna- 
tive presentations with essentially the same final equa- 
tions can be produced in which the relaxation effect is 
mainly concentrated in the substrate, or in both. For 
concreteness, we refer only to the previous realization. 



III. CONSTANT VELOCITY RESULTS 

In this section we solve numerically Eqs. ([3]) and ^ 
for models A and B, at constant values of the driving 
velocity and in the limit of 770 — >■ 0. This is achieved 
through the following protocol. We set the position of 
the driving at some constant value C/, and iterate the 
evolution equations using a simple first order algorithm, 
until an equilibrium position is obtained for x (i.e, until 
X is sufficiently close to zero). This process is assumed 
to occur at fixed time, as 770 — >■ 0. After this, time is 
advanced a small dt, which produces an increase of [/ — > 
\J -f Vdt, and the process is repeated. 

The disordered pinning potential we use is obtained 
through the following interpolation mechanism. For all 
points of the horizontal axis with coordinates of the 
form Xi = ai, with i an integer, a Gaussian random 
number 0^ of zero mean and unitary variance is gen- 
erated. Each produces a pinning potential that is 



zero outside the interval [a;i_2, 2:^+2], and has the form 
Wi{x) = hQ(t)i [cos (7r(a; — Xi)/2a) -\- 1] inside it. The to- 
tal W{x) is the sum of the individual contributions: 
W{x) = 'Yl,iWi{x). Note that there are at most four 
non zero terms in this sum. In the case of model B, we 
scale down the forces /i and /2 by a factor \/2, in order to 
have typically the same overall friction force as in model 
A. Also, we will present all results in dimensionless form, 
using a as the unit of length, rjo? /F^hQ as the unit of 
time, and Ff^h^la as the unit of force. 
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FIG. 3: Average friction force as a function of the pulling 
velocity, for the two models sketched in the previous figure. 
Both of them display velocity weakening (results here and 
in the following figures were obtained for fco = fci = 0.3). 
Model V also shows a velocity increase of friction force at low 
velocities that follows an a -I- (3V^^^ law (dotted line). 

The results for the average friction force as a function 
of velocity are contained in Fig. |3l The two models 
clearly display velocity weakening in some velocity range 
(model B also displays increase of friction force at the 
lowest velocities, we address this effect later on). The 
qualitative justification for the existence of VW applies to 
the two models in the same way, and goes as follows. At 
large velocities, the driving induces temporal variations 
in the dynamics with a period T = a/V . On the other 
hand, the dissipative element possess a time constant r 
given hy T = rjk (k = + ki for model A, and k — ki 
for model B). For t ^ T the ki spring is effectively 
blocked by the dash-pot, namely 9 is nearly constant for 
sufficiently large velocities. On the contrary, at the lowest 
velocities, there is plenty of time for 9 to adapt, which 
generates an additional dissipation that manifests in an 
increase of the friction force. 

The time evolution of the instantaneous friction force 
reveals some differences between models A and B. Let us 
start by describing the evolution for model A. The in- 
stantaneous friction force at different velocities is shown 
in Fig. 01 As previously stated, the evolution for large ve- 
locities corresponds to that in the original FT model with 
a driving spring /cq. We observe the typical jumps in the 
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FIG. 4: Instantaneous friction force as a function of time, 
for model A (see Fig. EJa)), at different pulling velocities. 
At large velocity (upper panel) the model behaves almost as 
a normal PT model with spring constant ko. As velocity 
is reduced, there is an additional response associated to the 
dissipative element rj and the ki spring. At the lowest ve- 
locities (lower panel), and except for time intervals of order 
r = {ko + ki)ri after the force jumps, the response is equivalent 
to a normal PT model with a pulling spring of effective con- 
stant fee// = fcofei/(feo + fei). In time intervals of order t after 
the jumps, ASs can be observed. One example is highlighted. 



friction force, corresponding to rapid rearrangement of 
the X variable, that is the origin of a finite average force. 
At the lowest velocities, we see that after the jumps, there 
is a rapid recovering of the force, which, with a time con- 
stant given hy T — (ko + ki)r] approaches a linear behav- 
ior in which force increases with a spring constant corre- 
sponding to an effective stiffness k^ff = koki/{kQ + ki). 
So there is a transition in the stick stage of the movement, 
passing from an effective constant feo large velocities, 
to keff at low velocities. 

In the case of model B (Fig. [5]), although velocity 
weakening is equally observed, the time evolution of in- 
stantaneous force does not display a transition between 
two different spring constants, instead, it always corre- 
sponds to an evolution with stiffness feg- This difference 
in the time evolution of friction force can be a potentially 
useful tool to discern between the different mechanisms 
producing VW in a particular experimental realization. 

In addition to display VW, both A and B models also 
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FIG. 5: Same as previous figure for model B (see Fig. [IJb)). 
Although this model also shows VW and ASs, the slope of 
the force increase during the 'stick' stage of the movement 
corresponds always to a stiffness feo. 

possess aftershocks, as we will discuss now. Let us con- 
centrate for concreteness in model A. The jumps in the 
evolution of the friction force correspond to transitions 
between different metastable states of the system. As it 
was already discussed, and it is clear in Fig. HI at very low 
velocities the instantaneous friction force after a jump 
displays a dip, and then a recovery during a typical time 
T. If this force increase is enough to destabilize the new 
position in which x is located, this will produce a second 
jump to a more stable minimum. This second jump is an 
aftershock to the first one. Note that aftershocks (that 
can in fact be more than one for a single initial jump) 
occur in a time scale of order r which depends exclu- 
sively on the internal dynamics of the system. In fact, at 
sufficiently low velocities, it is observed that sometimes 
there exist sequences of a few jumps that occur close in 
time (actually, within a time interval of order r) and they 
represent the initial shock and its aftershocks. Some par- 
ticular examples are highlighted in Figs. |4]and[5] A more 
quantitative description of ASs can be made by generat- 
ing the histogram of time intervals between consecutive 
jumps in the instantaneous friction force. The typical 
time interval expected between main jumps is given by 
^ a/V. On the other hand, the typical time interval 
for the appearance of an aftershock is r. It is clear that 
in order to unambiguously identify ASs we must have 
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V <C a/r. This is in fact what is observed in the results 
shown in Fig. IH^a). At large velocities, we see a smooth 
distribution of time intervals that for the present form 
of the pinning potential can be rather well approximated 
by a Gaussian distribution. In the absence of relaxation 
(i.e., rj oo), as velocity is reduced we would simply 
obtain the same distribution but stretched to larger time 
intervals. These rescaled curves are shown in Fig. [Blja) 
as continuous lines. We see that for sufficiently low ve- 
locities, the true distribution of time intervals displays a 
peak at small time intervals of order r (independent of 
V). This is the peak corresponding to ASs. For the pa- 
rameters used in Fig. [6Ka) we find that about one quarter 
of all jumps are ASs. We emphasize that ASs are visi- 
ble only at sufficiently low velocities, otherwise they are 
masked by the jumps that are produced by the driving 
itself. 
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FIG. 6: (a) Histograms of time intervals between jumps for 
model A at different driving velocities. Al large velocities, the 
distribution is well fitted by a Gaussian distribution (contin- 
uous line). When this distribution is rescaled to be applied 
to lower velocities (dotted lines) it only fits the large At part 
of the distributions. The additional events with low At are 
ASs (for the lowest velocity, ASs are indicated by the dashed 
region). Note that ASs are not affected by V, once V is small 
enough to allow their observation, (b) The same qualitative 
behavior is observed for model B, although with slightly dif- 
ferent forms of the distributions. 

The situation for model B is quite similar. The his- 
tograms of time intervals (shown in Fig. IHl^b)) display 
the same features, although the actual form of the distri- 
butions are slightly different. In the case of model B, ASs 
occur when in the process of relaxing the variable 0, the 
system finds an instability point for the center of mass 
variable x. From Fig. [5] we can conclude that the aver- 



age time between main shocks a/V has to be larger than 
about 10 in order than ASs can be identified. Restoring 
units, this gives V < O.lFpfho/rja as necessary to observe 
ASs. 

Model B presents at the lowest velocities a region of ve- 
locity strengthening (see Fig. namely , an increase of 
friction force upon velocity increase. A detailed analysis 
of the time evolution of x and 9 reveals the mechanism 
that originates the effect. Consider the limit V 0, 
and imagine the trajectory of the coordinates of the sys- 
tem in the x-9 plane. The trajectory consists of smooth 
pieces, and abrupt transitions when instabilities occur. 
At a somewhat larger value of V, the trajectory follows 
basically the same path that in the V limit, but 
the abrupt transitions between smooth pieces of trajec- 
tories is slightly delay due to the time constant of the 
9 variable. This produces an increase in friction force 
through the same mechanism that operates in Eq. ([!]) 
for non-zero rjQ. In that case an increase of friction force 
as ~ 1/2/3 jg ^qY[ knownl^'-"*, and originates in the delay 
that the x variable experiences in trying to follow the 
quasistatic behavior, when 77 > 0. In the present case, 
it is the 9 variable that takes some time to adapt to the 
evolving landscape. In this case however, the instability 
occurs through a Hopf bifurcation, and this modifies the 
velocity dependence, which now is of the form V^^^. In 
fact, the low velocity increase of the friction force is very 
well fitted by this velocity dependence (see Fig. 



IV. RESULTS FOR SLIDE-STOP-SLIDE 
EXPERIMENTS 

The mechanisms that lead to VW in geological mate- 
rials produce also characteristic features in experiments 
in which the velocity of the driving is changed in some 
abrupt way. For instance, if driving is stopped, there 
is typically a logarithmic decrease in time of the stress 
in the system, as the internal degrees of freedom relax. 
When driving is reinitiated, there is a transient in which 
the stress displays an overshoot with respect to the sta- 
tionary value, which reflects the higher energy barrier 
that has to be overcome to reinitiate sliding, after a pe- 
riod of rest. It is interesting to see if this phenomenology 
is also observed in the simple generalizations of the FT 
model we are discussing here. Due to the large varia- 
tions of the instantaneous friction force in our system, it 
is clear that a systematic effect will only be observable if 
an average over many realizations is performed. In Fig. [7] 
we show this average for the instantaneous friction force 
as a function of time in a system that is initially driven at 
a fixed velocity. The driving is stopped at time and re- 
initiated at time ti with the same original velocity. The 
results show that the stress reduction in time after sliding 
stops, and the peak after reinitiation, are qualitatively 
reproduced. The detailed form of the time dependencies 
however, cannot be expected to match the experimental 
one o^'^° , as they depend on collective effects that are not 



7 



captured by the simple models discussed here. 
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FIG. 7: Average friction force over many realizations for a 
system tiiat is initially sliding at velocity V. Sliding stops at 
t = t\ and is reinitiated at the same v &X t = t2- Results are 
for model A [(a) V = 0.1, 2000 realizations] and [(b), V = 100, 
10^ realizations] with other parameters as in previous figures. 



V. SUMMARY AND CONCLUSIONS 



phenomena may occur and be observable in a variety of 
situations. VW and ASs were obtained here by introduc- 
ing two different types of relaxation in the standard PT 
model. In one of them, the driving spring shows a vis- 
coelastic response. In the second, the contact region was 
assumed to have some spatial extent and to be described 
by more than a single variable. In both cases, there is an 
internal time constant that characterizes the response of 
the system. VW and ASs are observed when the driving 
velocity leaves enough time to the contact to adapt its 
internal degree of freedom as sliding proceeds. 

Temperature effects have not been consider in our anal- 
ysis. It is clear that thermal fluctuations will be a com- 
peting effect against the observation of VW ans ASs. As 
thermal fluctuation effect are more important in smaller 
size devices, the observation of ASs with the AFM may 
be a challenging task. 

The two model analyzed here were constructed having 
in mind situations of very small contact regions, like those 
realized in the surface force apparatus, or even at smaller 
scales, with the AFM. Applications to larger contact 
regions likely require the consideration of an extended 
model with many degrees of freedom. In this respect the 
consideration of a model similar to the present model B in 
a situation with many contact point has been considered 
for instance inii^, and in fact led to a phenomenology 
comparable with that observed in earthquakes. The con- 
sideration of model A in a similar context has not been 
done yet, and remains as a future prospect. 



Velocity weakening and aftershocks are two key ingre- 
dients that characterize the dynamics of seismic phenom- 
ena. They are deeply related with the existence of re- 
laxation effects within the sliding materials. When going 
progressively to smaller spatial scales, these two phenom- 
ena have been difficult to observe. 

We have shown however that some simple relaxation 
mechanisms, that may reasonable be acting in microcon- 
tact experiments, can lead to the expectation that these 
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